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We present results on the connection between the vorticity equation and the shape of the single- 
point vorticity PDF. The statistical framework for these observations is cast in form of conditional 
averages. The numerical evaluation of these conditional averages provides insights into the intimate 
relation of dynamical effects like vortex stretching and vorticity diffusion and non-Gaussian vorticity 
statistics. 
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One of the main goals of theoretical turbulence re- 
search is to develop a statistical non-equilibrium theory of 
turbulence from first principles. While this goal has not 
been achieved so far, many phcnomcnological theories ex- 
ist, explaining a variety of phenomena of turbulent flows 
P, For example, the classical K41 [3] theory man- 
ages to predict the observed energy spectrum, and later 
theories like K62 or the multifractal model [H even give 
arguments for the intermittent scaling of the structure 
functions for the velocity increments. However, all these 
classical theories do not directly link to the equation of 
motion of a fluid, i.e. the Navier-Stokes (or equivalently 
the vorticity) equation. 

There are, however, works attempting to derive a di- 
rect statistical theory of turbulence from first principles. 
For example, Lundgren Q derived an evolution equation 
for the turbulent velocity with the help of PDF methods. 
Similar results were obtained by Novikov |5|, however 
focusing on the vorticity. Extensions to turbulent com- 
bustion and many other applications were proposed by 
Pope Q. Mathematically all these approaches have to 
face the famous closure problem of turbulence, which is 
inherent to the highly nonlinear character of the equa- 
tions of motion. In Lundgren's work this problem results 
in an infinite hierarchy of evolution equations, involving 
an ever increasing number of velocities at different spa- 
tial points. In the works by Novikov a formal closure is 
achieved by introducing conditional averages, which en- 
ter the statistical equations as unknown functions. 

It is mainly due to experimental efforts and direct nu- 
merical simulations of turbulence that today we know 
that turbulent flows are governed by coherent structures. 
In fully developed, homogeneous and isotropic turbu- 
lence these structures appear as filamentary vortex tubes. 
Each of these tubes generates a swirling velocity field. 
The ensemble of vortices present in a turbulent fluid then 
interacts in a nonlinear manner, leading to the complex 
spatio-temporal structure of turbulent flows. In the inter- 
action of these structures, some basic dynamical mecha- 
nisms can be identified: advection, vortex stretching and 
viscous diffusion of vorticity. 

It is a scope of the present work to establish a link be- 
tween these dynamical aspects of turbulence and the non- 



Gaussian vorticity probability functions (PDFs). To this 
end, we present an evolution equation for the turbulent 
vorticity PDF in the spirit of Lundgren and Novikov. The 
appearing conditional averages can be evaluated with the 
help of highly resolved direct numerical simulations of the 
vorticity equation. The results help to reveal the connec- 
tions between turbulent dynamics, coherent structures 
and non-Gaussian statistics. 

The temporal evolution of the vorticity i^{x, t) is gov- 
erned by the vorticity equation, 

— +u -S/u: = S -uj + vAu} + F, (1) 

where u{x,t) denotes the velocity field and Sij = 
i (^^{x,t) + ^^{x,t)j denotes the rate-of-strain ten- 
sor. As we want to focus on incompressible fluids, the 
velocity field can be obtained from the vorticity field via 
Biot-Savart's law. v denotes the kinematic viscosity and 
F{x,t) is an external forcing in order to achieve a sta- 
tistically stationary state. Starting from the fine-grained 
PDF f{n;x,t) = 6{u:{x,t) ~ standard PDF methods 
yield a kinetic equation for the turbulent vorticity PDF 

fin;x,t) = {fin;x,t)) iii, 

^/ + V,-{(m|J^)/} = -Vo-{(S'-u; + i.A,a; + F|n)/}. 

(2) 

The right hand side of this kinetic equation for the vor- 
ticity PDF reveals the different dynamical influences: 
the average vortex stretching term, vorticity diffusion 
and the forcing conditioned on the sample-space vortic- 
ity. Taking into account homogeneity, the advective term 
vanishes, as both the conditional average as well as the 
PDF do not depend on the spatial coordinate. It is ar- 
gued in [s] that the conditional average of the forcing 
decays with increasing Reynolds number like {F\Q,) ~ 
Re~^^^. Given a sufficiently high Reynolds number, the 
large-scale forcing should not affect the smallest scales of 
the fiow. As the coherent structures live on these scales, 
a negligible infiuence of the forcing is physically sound. 
With these simplifications the kinetic equation reads 

^J = -Wn-{{S-u: + ,^AMf^)f}. (3) 
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In order to achieve a statistically stationary state, the 
right hand side of this equation has to vanish. This is 
only possible by a statistical cancellation of the appear- 
ing terms. In the following we will, in addition to ho- 
mogeneity, consider isotropic turbulence. This imposes 
further constraints on the statistical quantities. It fol- 
lows that {iyAu}\n) - n. With {S ■ u}\n) = {s\n) ■ rj, 
the conditioned rate-of-strain tensor can accordingly be 
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written down as {Sij\fl) = g{n)6ij — 3h{n)-^ 
two scalar functions h and g only depending on the ab- 
solute value of the sample-space vorticity. The trace of 
this conditionally averaged tensor has to vanish, which 
yields g = h. With this it is easy to show that fl is an 
eigenvector of {Sij\fl), {Sij\fl)Uj = —2g{fl)fli with the 
eigenvalue Ai ~ —2g{Q). The remaining eigenvalues can 
directly be determined due to the trace condition and 
isotropy, A2,3 = g{^)- Thus a cancellation of the two 
terms on the right hand side of ([3]) is possible in a statisti- 
cal sense. While these observations show, which dynami- 
cal effects statistically have to cancel, these equations do 
not suffice to directly calculate the vorticity PDF. This 
can be achieved by taking into account the homogeneity 
of the flow. Calculating the Laplacian of the vorticity 
PDF yields 
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(4) 

In case of a single component, this equation reduces to 
= "MT^^^^I^a;)/ - ■^{{Vu)x)'^\VL^)f, which suffices 
to determine the functional form of the vorticity PDF for 
homogeneous (and not necessarily stationary) flows 0. 

An appealing kinetic equation arises, when combining 
the kinetic equation with the homogeneity relation 
For homogeneous turbulent flows, the temporal evolution 
of the vorticity PDF can then be described by 
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For a single component, equation ([5]) reads 
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which yields the stationary solution 



{u{VLO,f\n^)f, (6) 
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with a normalization constant TV. This shows, that in a 
stationary homogeneous flow the vorticity PDF is deter- 
mined by the dynamical effect of vortex stretching and 
the vorticity gradient. 

Further information on the shape of the PDF is ob- 
tained by studying the non-stationary case of equa- 
tion ([3]) with the method of characteristics. The ordi- 
nary differential equation for the characteristics reads 
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512 164 0.082 0.0001 2.19 26.73 0.008 0.63 1.6 



Table I: Major simulation parameters. Number of collocation 
points TV, Reynolds number based on the Taylor micro-scale 
R\, root-mean-square velocity Wrms, kinematic viscosity ly, in- 
tegral length scale L, large-eddy turnover time T, Kolmogorov 
length scale rj, Kolmogorov time scale r^, kmaicil characterizes 



with the spatial resolution of the smallest scales. 



tt = {S ■ uj + vAxUjirt) with the solution describing the 
temporal evolution of the sample-space vorticity. 

In the following we will study the PDF /(r2^;t) of a 
single component of the vorticity. Equation ^ then 
reduces to 
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The statistical evolution of the sample-space vorticity fl^ 
is accordingly given by = {(S ■ lj)x + J^AajCJa;!!!^;), i.e. 
by the sum of the conditionally averaged vortex stretch- 
ing term and the conditionally averaged vorticity diffu- 
sion. These averages will be evaluated numerically below. 

The turbulent fields under consideration in the present 
work are generated by a standard, dealiased Fourier- 
pseudo-spectral code [8|,[§| for the vorticity equation. The 
integration domain is a triply periodic box of box-length 
2TT. To obtain a statistically stationary flow we apply a 
large scale forcing. The time-stepping scheme is achieved 
by a third order Runge-Kutta scheme ■ 

For the present work, we conduct two different types 
of simulations. For an estimation of the conditional aver- 
ages determining the stationary PDF a run in the statisti- 
cally stationary regime is performed, table |T] sums up the 
major simulation parameters. In order to gain deeper 
insights into the formation of the PDF a second, non- 
stationary run with comparable simulation parameters 
is performed. The initial condition of this run exhibits 
a Gaussian vorticity distribution and the same energy 
spectrum as the statistically stationary simulation. The 
held evolves under the dynamics of the vorticity equation 
and during the relaxation to the statistically stationary 
state the non-Gaussian vorticity PDF emerges. 

Regarding the stationary situation, figure [1] shows the 
numerically evaluated conditional averages of equation 
([8|). The conditionally averaged vortex stretching term 
is positively correlated with the vorticity component, 
while the conditionally averaged vorticity diffusion shows 
strong anticorrelations. These tendencies can be physi- 
cally understood; in presence of strong vorticity the vor- 
tex stretching term causes a self-amplification of vortic- 
ity, the diffusive term then tends to deplete this vortic- 
ity. Checking the validity of equation ([7]) with our 
numerical data, the comparison of the directly estimated 
and reconstructed vorticity PDF is depicted in figured) 
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Figure 1: Conditional averages of the vortex stretching term 
a) ((S ■ Lo)x\^x), the vorticity diffusion b) {v1^lox\^^) and 
the sum of both terms. The sum vanishes, as required for 
statistical stationarity. The squared vorticity gradient c) 
{v{y ujx^'^\^x) exhibits a strong dependence on f2a;. 




Figure 2: Logarithmic plot of the vorticity PDF estimated 
directly from our data and the reconstructed PDF according 
to eq. ((6}. The agreement over several orders of magnitude 
is excellent, slight deviations are visible in the far tails of the 
PDF. 



An excellent agreement over seven orders of magnitude 
is given. As equation ([6]) yields a Gaussian in case 
of {(S ■ ijS)x\^x) and {i'{'Vujx)'^\^x) ~ const., 

the highly non-Gaussian shape of the vorticity PDF 
can be tracked down to the strong fJjj-dependence of 
{{iy\/ujx)'^\^x)- Further motivation for the functional 
form of the conditional averages can be given when ex- 
amining the basic structures present in the flow. As the 
visualizations in Fig. [3] suggest, the flow consists of elon- 
gated vortex tubes, which can be modeled with Burgers 
A Burgers vortex exhibits a vorticity field 



vortices 
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according to u; = uj{r)ez = characterized 
by the strain parameter a and circulation F. For this 
structure the vortex stretching term and the squared vor- 
ticity gradient are readily calculated to • a; = au) and 




Figure 3: Volume rendering of the absolute value of vorticity 
above a fixed threshold for different stages of the simulation 
(from top left to bottom right: initial condition, 0.11 T (200 
timesteps), 0.38 T (600 timesteps), 3.53 T (6000 timesteps)). 



(Vw^)2 = {^Lj) = ^u}^- That means, for fixed a, F 
and the vortex stretching term is a linear function of 
the vorticity, whereas the squared vorticity gradient turns 
out to be a quadratic function of the vorticity. Think- 
ing of turbulence as an ensemble of Burgers-like vortices, 
this picture already captures the main features of the 
conditional averages shown in Fig. [1] Deviations from 
this simple argument are possible due to the fact that 
the circulation of a vortex tube in a turbulent flow is 
not independent of the surrounding strain field and that 
straight vortex tubes are not the only structures present. 

While helping to explain the non-Gaussian nature of 
the turbulent vorticity, the above considerations give no 
direct information on how the dynamical processes of vor- 
tex stretching and vorticity diffusion go along with the 
formation of the characteristic shape of the PDF. To elu- 
cidate this issue, we turn to the non-stationary simula- 
tions starting with a vorticity field exhibiting a Gaussian 
distribution. Under the temporal evolution of the vor- 
ticity equation ^ strong spatio-temporal correlations in 
form of e.g. vortex tubes are generated and the initial 
Gaussian distribution relaxes to the non-Gaussian PDF 
observed for the statistically stationary regime. The tem- 
poral evolution of the vorticity field is visualized in fig- 
ure [H The initial condition with the Gaussian PDF ap- 
pears unstructured, but already at 0.11 T the emergence 
of small vorticity worms can be observed. In the course 
of the simulation these structures grow stronger, after 
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Figure 4: Illustration of the deformation of the PDF for 
the non-stationary simulation. Upper figure: The vortic- 
ity PDF for O.llT is already non-Gaussian, yet has not 
developed strong tails. Lower figures: conditional sum 
{{S ■ (jl})x + fAaiLLixlQx) shown for two different y-ranges. The 
inner region of the vorticity PDF is quenched together due 
to dominant vorticity diffusion while the outer regions are 
stretched towards higher vorticity values due to dominant vor- 
tex stretching. 



0.38 T thin vortex tubes can be observed. Interestingly 
these tubes tend to cluster, as can be seen in the snapshot 
taken at 3.53 T. 

The temporal evolution of the sum of the conditional 
averages of eq. dS]) is evaluated for the non-stationary 
situation. An example for O.llT is depicted in the lower 
part of figure ID As expected, the sum of both aver- 
ages does not cancel like in the statistically stationary 
run. While the initial condition is purely diffusion domi- 
nated as the vortex stretching has not started generating 
structures, the formation of two distinct regions can be 
observed in the course of the simulation. Low absolute 
values of the vorticity are dominated by the conditional 
vorticity diffusion while larger absolute values are ampli- 
fied due to a dominant vortex stretching. These differ- 
ent regions remain as the PDF is approaching station- 
arity, but the two terms tend to cancel more and more. 
After reaching the statistically stationary state the sta- 
tistical balance of the conditional averages is recovered. 
How the unbalanced conditional sum affects the forma- 
tion of the pdf may be interpreted with the method of 
characteristics from figure[31 The inner, diffusively domi- 



nated region is quenched towards zero vorticity, while the 
outer, vortex stretching dominated regions are stretched 
towards larger absolute values of vorticity. The emerging 
physical picture fits well to the common understanding 
of turbulent flows. The stretching of the PDF towards 
larger absolute values of vorticity corresponds to the self- 
amplification of vortex filaments due to vortex stretching. 
The lower-valued vorticity, which corresponds to more 
unstructured regions of the flow is depleted by the diffu- 
sive term. Hence the formation of non-Gaussian vorticity 
PDF can be tracked down to the two physical mecha- 
nisms of vortex stretching and vorticity diffusion. 

To summarize, we reported on theoretical and numer- 
ical results on the link between the vorticity equation, 
coherent structures and the non-Gaussian distribution 
of vorticity. Based on classical works by Lundgren and 
Novikov, an investigation of the kinetic equation of the 
one-point vorticity PDF reveals that the conditional aver- 
ages of vortex stretching and vorticity diffusion determine 
the temporal evolution and shape of the vorticity PDF. 
A closed expression for the stationary vorticity PDF 
was found in terms of the conditional averages of vor- 
tex stretching and squared vorticity gradients. Numeri- 
cal simulations confirm this relation with a high degree 
of precision. Further investigations of a non-stationary 
flow reveal that during the transition to the stationary 
state, two distinct regions of the vorticity PDF can be 
found: the inner region of this PDF is quenched due to 
dominant vorticity diffusion while the development of the 
fat tails can be associated with the stretching of strong 
vortices. 

Hence this work highlights a direct connection between 
basic dynamical features of turbulence and their statisti- 
cal consequences. These results encourage to study tur- 
bulent flows in terms of coherent structures as a main 
source for non-Gaussian statistics. 

We acknowledge insightful discussions with F. Jenko, 
O. Kamps, M. Vofikuhle and A. Daitche. Computa- 
tional resources are allocated at the LRZ Munich (project 
h0963) and on the BOB cluster at RZ Garching. Volume 
rendering produced with VAPOR, www . vapor . ucar . edu. 



[4] 
[5] 
[6] 



Electronic address: l mwilczek@uni-muenster.del 
U. Frisch, Turbulence - The Legacy of A.N. Kol- 
mogorov (Cambridge University Press, Cambridge, Eng- 
land, 1995). 

R. Friedrich and J. Peinke, Physical Review Letters 78, 
863 (1997). 

A. N. Kolmogorov, Proceedings of the Royal Society: 
Mathematical and Physical Sciences (1990-1995) 434, 9 
(1991). 

T. S. Lundgren, Physics of Fluids 10, 969 (1967). 
E. A. Novikov, Fluid Dynamics Research 12, 107 (1993). 
S. Pope, Turbulent Flows (Cambridge University Press, 
Cambridge, England, 2000). 



5 



[7] E. S. C. Ching, Phys. Rev. E 53, 5899 (1996). [10] C. Shu and S. Oshcr, ,J. Comp. Phys. 77, 379 (1988). 

[8] C. Canuto, M. Hussaini, A. Quarteroni, and T. Zang, [11] J. Burgers, Advances in applied mechanics, Academic, 

Spectral Methods in Fluid Dynamics ( Springer- Verlag, New York (1948). 

Berlin, 1987). 

[9] T. Y. Hou and R. Li, J. Comp. Phys 226, 379 (2007). 



